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UNM  Raport  No.  ME- 145 (68)  AF08R-993-3 


1.0  OBJECTIVES 


The  prtMnt  research  aims  toward  developing  structural  theo¬ 
ries  which  can  be  used  to  diagnose  the  fracture  damage  of  struc¬ 
tures  and  to  assess  the  reliability  of  the  damaged  structures. 
Such  theories,  once  verified  expert sent al ly,  will  be  available 
for  design  engineers  to  apply  to  the  damageable  structures.  It 
can  also  be  the  base  for  Methodology  which  analytical  and  test¬ 
ing  engineers  can  develop  to  diagnose  and  assess  the  reliability 
of  the  existing  structures.  It  is  assumed  that  the  structures 
under  consideration  aay  develop  damage  through  extreme  excita¬ 
tions.  Such  damage  can  be  defined  as  cracks  occurred  in  the 
structures,  the  amount  of  energy  dissipation,  the  deformation  or 
any  combination  of  the  three.  Therefore,  it  is  important  to 
know  when  damage  has  occurred  in  a  structure.  When  it  has,  it 
is  desirable  to  be  able  to  locate  it  and  estimate  its  extent. 

Following  the  damage  diagnosis,  the  principal  components 
for  damage  assessment  are  (1)  to  develop  theories  for  the  relia¬ 
bility  assessment  of  damageable  structures  and  estimate  the  dam¬ 
age  in  a  structure,  and  <2)  to  develop  and  improve  mathematical 
models  which  simulate  the  behavior  of  damageable  structures. 

These  assessment  theories  will  assist  the  engineers  to  achieve  a 
specified  predictive  accuracy  in  design,  and  to  obtain  a  more 
realistic  assessment  of  the  existing  damageable  structures. 


2.0  RESEARCH  PLANNINBi 


The  rttwrch  effort,  sponsored  by  the  Air  Fore*  of  Scienti¬ 
fic  Research,  eade  progress  during  the  research  years  1981-1984 
to ward  developing  a  consistent  structural  theory  in  daeage  diag¬ 
nosis  and  reliability  asseeseent  Cl  through  163  (reference  to 
Section  10.0,  publications  by  the  principal  investigator  F.  D. 
Ju>.  The  accompli sheen ts  are  in  phenomenal  correlation,  in 
structural  modeling  and  in  quantitative  diagnosis  of  fracture 
damage.  Deterministic  analysis  mas  the  goal  if  possible.  In 
real  life,  the  excitation  and  even  the  structural  characteris¬ 
tics  can  be  random.  Moreover,  the  measurements  may  be  inadequa¬ 
te,  inaccurate  or  garbled  with  noise.  The  reliability  assess¬ 
ment  will  have  to  be  probabilistic.  With  the  establishment  of  a 
practical  theory  in  mind,  a  plan  for  three  years  research  was 
proposed  in  1984  and  approved  beginning  January  1,  198S. 

2.1  icaaiffitf  j  bU  i  t  v  ..Thtar.  y » 

The  transmissibility  is  defined  as  the  ratio  between  the 
response  and  the  excitation  at  two  different  stations  on  a  struc¬ 
ture.  Cracks  soften  the  structure  and  thus  cause  a  change  in 
transmissibility.  The  locations  of  transducers  which  measure 
the  transmissibi 1 ities  as  well  as  the  locations  of  excitation 
have  been  optimized  for  the  necessary  number  of  transducers  and 
for  the  maximum  measurement  in  the  change  of  transmissibi  lity. 
Damping  property  in  the  structure  is  included  in  the  model.  The 
research  would  propose  an  optimal  rule  to  locate  the  response 
transducers  and  the  number  of  measurements  necessary  in  isola¬ 
ting  the  location  of  a  crack  damage  and  in  assessing  the  inten¬ 
sity  of  the  damage. 
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The  analytical  researches  on  the  modal  methods  o-f  damage 
diagonal  a  <hnfh  the  modal  -frequency  theory  and  the  transnussibi- 
lity  theory)  are  based  on  the  theoretical  modeling  of  the  crack 
with  a  fracture  hinge  by  the  principal  investigator  Cl ,2, 3, 8, 9, 
11].  No  experiment  was  then  available  to  establish  and  to  veri¬ 
fy  its  existence.  This  phase  of  research  was  considered  moet 
critical  in  the  decision  to  continue  the  research.  For  the  com¬ 
pletion  of  tho  experimental  reeaarch,  the  additional  equipment 
and  material  supports  -from  the  Mechanical  Engineering  Depart¬ 
ment  t  UNM  and  the  Sandia  National  Laboratories  are  gratefully 
acknowl edged . 


The  diagnostic  theories  are  essentially  inverse  problems, 
for  which  the  solution  say  not  exist  or  aultiple-vaiued.  In  the 
direct  problee,  the  structure  and  the  dsesge,  the  location  and 
the  extent  of  fracture,  are  given.  The  pre-  and  the  post-damage 
structure  can  be  analyzed  for  its  corresponding  frequency  equa¬ 
tions.  The  pre-  and  post-damage  dynaeic  characteristics  are, 
therefore,  enumerable.  However,  in  the  inverse  problee,  the  dam¬ 
age  is  not  known  a  priori.  The  sufficiency  of  measurements  is 
always  an  uncertainty.  Furthereore,  as  the  damage-softened 
structure  tends  to  lower  the  eodal  frequencies,  the  neighboring 
modal  frequencies  could  very  well  experience  s  cross-over .  In 
the  modal  frequency  theory,  such  cross-over  cannot  be  ascer¬ 
tained.  The  types  of  uncertainties  have  to  be  address  in  order 
to  caution  the  users  of  the  theory  and  to  suggest  naw  theories 
to  improve  the  accuracy  in  using  ths  theory. 
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.4  Probabilistic  Theory  of  Multiple  Fracture  Pisiributipnt 

The  dtMB*  diagnosis  for  a  singla  crack  has  been  developed 
succsss-f ul ly .  Ho— vt ,  givsn  an  external  randoa  excitation,  the 
cracking  location  cartainly  is  not  a  fixed  pattern.  Tha  random¬ 
ness  of  the  first  cracking  location  eay  be  correlated  to  the  mat¬ 
erial  of  the  e  structures  and  the  characteristics  of  the 
external  excitation.  If  the  excitation,  such  as  earthquake,  is 
extre— ,  it  is  anticipated  that  multiple  cracks  mill  occur.  The 
occurrence  and  location  of  the  second  crack  may  be  correlated  to 
the  first  crack,  and  so  on.  In  the  modal  diagnostic  theory, 
because  the  a  priori  knowledge  in  the  number  of  cracks  is 
lacking,  the  diagnosis  is  uncertain.  Accordingly,  it  is 
necessary  to  develop  the  probability  distribution  of  the 
cracking  location.  The  established  probability  distribution  of 
the  cracking  locations  can  help  to  identify  multiple  cracks  at 
least  in  the  probabilistic  sense. 


2.5  Modeling  of  Elastic-Plastic  HultlrfleqrTPf.-frfdBtt 

The  research  is  the  completion  of  the  study  beginning  with 
the  Single-Degree -of -Freedom  <8DF>  system  14,9,7,10,14,161.  The 
result  establishes  a  model  for  a  generic  class  of  damaged  struc¬ 
tures.  With  such  a  nonlinear  MDF  stochastic  model,  the  response 
of  a  damaged  structure  subjected  to  arbitrary  random  excitation 
can  then  be  analyzed.  The  procedure  of  characterizing  the  sys¬ 
tem  and  its  governing  stochastic  differential  equation  results 
in  uncoupled  differential  equations  of  its  mean  response  and  its 
random  component.  The  responses,  their  moments  and  crossmoments 
are  solved  for  the  statistical  — asures  of  the  cumulative  defor¬ 
mation  and  the  cumulative  energy  of  dissipation  as  damage  assess¬ 
ment. 


2.6  Reliability  of  Structures  with  Stiffness  and  Strength 

Ptaratiitipni 

Whan  a  structure  of  friction  materials,  such  as  concrete, 
is  subjected  to  strong  randoe  excitation,  the  structure  eay 
undergo  inelastic  deforeations  during  certain  cycles  of  loading 
with  associated  cyclic  degradation  in  stiffness  and  strength. 

The  prediction  of  reliability  and  assesseent  of  daeage  depend 
upon  the  proper  eodeling  of  such  structures,  taking  into  consi¬ 
deration  of  the  randoe  characteristics  of  the  eaterial  variables 
as  well  as  the  excitation.  It  has  been  observed  experieental ly 
in  the  structural  degradation  eodels  that  the  rate  of  degrada¬ 
tion  eay  be  related  to  the  energy  dissipation  through  the  degra¬ 
ding  restoring  hysteretic  loop.  The  randoaness  of  the  stiffness 
and  strength  degradation  is  estieated  with  the  realization  of 
the  energy  dissipation.  The  reliability,  which  is  defined  by 
the  eaxieue  deforest ion  that  never  reaches  an  ass used  and  dete« — 
einistic  critical  level,  can  then  be  evaluated  based  on  the  es¬ 
tablished  theory.  However,  the  statistical  relation  between  the 
randoe  stiffness  degradation  and  the  energy  dissipation  is  an 
assueed  probability  density  function  since  currently  there  is  no 
data  available. 


-5- 


3.0  TRANSMI88IBIL1TY  THEORY* 


The  transaissibi X i ty  theory  use*  the  dynuic  characteris¬ 
tics  of  the  structural  responses.  Sines  the  excitation  aay  be 
carried  out  at  soee  specific  excitation  frequencies  or  so ae 
known  statistics  of  randoe  excitations,  the  uncertainties  aris¬ 
ing  through  either  insufficient  aeasurseents  of  aodal  frequen¬ 
cies  or  frequency  cross-overs  can  thus  be  avoided.  The  theory 
is,  therefore,  eost  suitable  for  fracture  diagnosis  in  coaplex 
structures.  The  analytical  eodel  has  the  fracture  characteris¬ 
tics  (e,4)  as  paraaetric  variables  as  in  the  aodal  frequency 
theory  (18, 19,303.  Hence,  the  general  theory  of  circuit  analogy 
(2,3,151  is  used.  The  vector  of  aodal  displaceaents  can  be  de- 
terained  froa  given  excitation  Fit)  •  exp (lot)  at  a  given 

excitation  station  (153  as, 

y  -  C  F  e.  C3. 1) 

With  the  known  set  of  shape  functions,  the  displaceaent  at  any 
point  in  the  pre-  or  post-daaage  structure  can  be  coaputed. 
Appropriate  locations  of  response  transducers,  placed  on  the 
structure,  will  yield  the  record  of  transaissibility  for  frac¬ 
ture  diagnosis. 


3.1  Transaissibility  Heaaureoenit 

Transaissibility  <T>  is  defined  as  the  aagnitude  of  the 
ratio  of  acceleration  at  a  response  station  on  the  structure  to 
the  force  applied  at  the  excitation  station  on  the  structure; 
that  is 


T  -  |w/F | ,  <3.2> 

where  the  displaceaent  *w*  aay  be  coaplex  to  take  into  considera¬ 
tion  of  structural  or  viscous  daaping.  To  facilitate  diagnosis. 


the  relative  treneeiesibi lily  tR^)  is  introduced  as  follows 

T  -  T  T 

Rt  -  —E— ^  -  1,  (3.3) 

o  o 

where  the  subscripts  c  and  o  indicate  with  and  without  crack 
respectively.  When  the  structure  is  excited  with  the  excitation 
(aeplitude,  frequency  and  location)  before  and  after  a  strong 
daeage— causing  excitation,  the  relative  transeissibi  lity  becoms 

R_  -  (y  /y  >  -  1  <3.4) 

I  CO 

where  y£  and  y^  are  the  aeplitudes  at  the  saee  response  station 
before  and  after  the  daeage.  The  relative  transeissibi lity  thus 
ranges  froe  -1  to  ®.  The  lieits  occur  at  the  post-  and  pre— dae¬ 
age  nodal  points  respectively  for  an  undaeped  structure.  Since 
structural  daeping,  however  weak,  always  exists,  it  is  possible 
only  to  identify  at  a  certain  frequency  a  point  of  seal  lest  aep¬ 
litude  of  oscillation.  A  pseudo -node-point  (PNP)  thus  defines  a 
point  of  locally  ainieua  aeplitude  of  lateral  oscillation  for  a 
specific  frequency.  The  relative  transeissibi 1 ity  is  a  function 
of  the  excitation  location  and  frequency,  the  location  of  the 
response  transducer  and  the  fracture  character i sties.  The  opti- 
eal  locations  of  excitation  and  response  are  decided  by  conve¬ 
nience  and  eaxieue  value  of  the  relative  transeissibi lity  for  as 
seall  a  crack  as  possible. 


3*2  Bl&aflQMt  SUUflfl1 

The  nueber  of  structural  eeebers  to  be  instrueented  for  res¬ 
ponse  transducer  is  governed  by  the  influence  range,  which  is  de¬ 
fined  by  the  proxieity  of  a  crack  to  a  response  station.  A  crack 
at  a  significant  distance  froe  a  transducer  is  said  to  be  out¬ 
side  of  the  influence  range  of  the  response  station.  In  a  fraee 
structure,  the  fraee  cell.  Joined  directly  to  an  instrueented 
eeeber,  is  called  a  neighboring  cell.  Any  crack,  developed  in  a 
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neighboring  cell*  Mill  effect  signif icently  the  relative  trans- 
mi  ssibility  of  the  instrumented  member.  The  influence  range, 
therefore,  should  include  only  the  neighboring  cells.  The  mini¬ 
mum  number  of  transducers  for  a  frame  structure  is  accordingly 


N  «  nm/2 


for  even  nm 


(nm+l)/2  for  odd  nm. 


(3.5) 


Hodal  shape  dictates  the  appropriate  choice  of  transducer  mem¬ 
bers.  For  instance,  at  the  fundamental  frequency  excitation, 
there  is  no  difference  in  choosing  any  mall  (or  column)  for  the 
response  station.  Hence,  Fig.  3.1  illustrates  a  three-story 
four — span  frame  Mith  possible  locations  of  response  transducers. 
Figure  2  shows  the  responses  measured  by  stations  (T)  on  members 
(16,21,22,24,27)  respecti vely ,  for  a  crack  on  member  17.  The 
curves  in  Figure  3.2  also  illustrates  the  insensitive  range  of 
the  crack  position  on  member  17. 


2S  26 


21 _ 22 

s 
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Fig.  3. 1  A  3-story ,4-span 
frame,  *T‘  denotes 
response  station 


Fig.  3.2  Relative  transmissi- 
bility  vs.  crack  loca¬ 
tion  (excitation  on  13) 


For  maximum  measures  of  relative  transmissibility,  it  is 
seen  that  the  response  station  should  be  near  the  pseudo-node¬ 
point.  However,  measurability  must  also  be  considered.  In 


order  to  avoid  misreading  the  measurement,  it  i»  recommended 
that  a  pair  O'f  offset  transducers  be  placed  at  each  response 
stat  ion. 


3.3  Damage  Diagnosis: 

The  relative  ‘frequency  eeasures  at  various  response  sta¬ 
tions  can  be  used  through  their  influence  ranges  for  the  loca¬ 
tion  of  the  cracked  eeeber .  For  instance,  when  the  response 
station  on  eeeber  21  records  a  relative  transeissibi 1 i ty  of 
2.127.  At  such  a  high  value,  it  can  rest  assured  that  the  other 
response  stations  will  register  relative  transeissibi 1 ity  less 
than  1.  The  cracked  eeeber  is  No.  21.  However,  if  the  value  is 
0.930,  the  cracked  eeebers  could  be  No.  5,  8,  20,  or  21.  Member 
S  could  be  ruled  out  since  softening  of  that  eeeber  would  left- 
shift  the  PNP  to  result  in  eost  likely  negative  values  for  the 
relative  transeissibi 1 ity.  For  other  eeebers,  the  records  of 
other  response  stations  can  help  to  identify  the  daeage.  For 
instance,  an  equal  reading  of  the  relative  transeissibility  on 
the  response  station  of  No.  22  ieplies  syeeetry.  The  cracked 
eeeber  is  then  the  coluen  No.  8.  If  the  relative  transeissi¬ 
bility  on  No.  22  response  station  is  small,  the  crack  will  be  on 
the  beae  No.  20.  Yet  an  increase  of  the  reading  on  No.  22  res¬ 
ponse  station,  the  damaged  eeeber  will  again  be  the  beae  No.  21 
itself  with  the  crack  at  a  less  sensitive  location.  There  is  no 
detereinistic  location  of  the  cracked  position  possible.  How¬ 
ever,  the  diagnostic  result  could  bring  the  region  of  daeage 
small  enough  to  use  some  local  diagnostic  eethod.  It  was  also 
shown  in  Fig.  3.2  that  there  is  a  range  of  crack  locations,  in 
which  the  crack  daeage  does  not  produce  sufficient  dynamic 
change  at  a  specific  excitation  frequency  to  allow  efficient 
diagnosis.  The  effective  ranges  of  a  3x4  frame  structure  are 
listed  in  Table  3.1  for  the  fundamental  frequency  excitation. 
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TABLE  3.2-Additional  Diagnos¬ 
able  Regions  at  the  Third 
Modal  Frequency 
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The  "hole"  can  be  filled  with  sow  higher  frequency  excitation. 
The  effective  ranges  of  the  ease  fraee  from  the  third  mode  fre¬ 
quency  excitation  are  listed  in  Table  2.  Of  course,  the  pseudo¬ 
node-points  of  those  instruwnted  sewbers  for  the  third  eode  are 
different  from  those  of  the  fundaeental  eode .  However ,  a  siailar 
procedure  say  be  taken  for  the  placeeent  of  the  corresponding 
transducers. 

The  transsissibility  theory  for  structural  fracture  diagno¬ 
sis  avoids  the  eultiple  frequency  problea  that  occurs  for  both 
of  the  last  two  uncertainties  in  the  eodal  frequency  theory  to 
be  explained  in  the  uncertainty  theory.  Art.  5.0.  The  theory 
proposes  the  use  of  a  single  excitation  frequency  at  both  pre- 
and  post— damage  measurement  of  responses.  The  wthod  is  limited 
only  by  the  number  of  transducers  to  be  placed  on  the  structure. 
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With  the  knowledge  of  the  resp onse-station  influence  range,  the 
damaged  member  can  be  mingled  exit.  It  is,  however,  to  be  pointed 
out  that  the  transaissibility  theory  does  not  in  itself  diagnose 
the  intensity  of  the  fracture  damage.  Some  combination  usage  of 
local  inspection  technique  or  even  the  modal  frequency  theory 
could  be  employed  for  local  and  intensity  identification. 
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4.0  EXPERIMENTAL  DIAGNOSIS  OF  FRACTURE  DAMAGE I 


The  investigation  aims  to  verify  experimentally  that  the 
structural  effect  of  a  cracked  section  can  be  represented  by  an 
equivalent  spring-loaded  hinge,  and  to  deaonstrate  experimen¬ 
tally  the  principle  of  fracture  diagnosis  in  structures  by  the 
modal  frequency  theory  £1,113. 


4. 1  The  Modal  Frequency  Theory* 


The  modal  frequency  theory  uses  the  structural  dynamic  con¬ 
cept  that  the  modal  frequencies  change  as  a  result  of  structural 
softening  from  the  presence  of  a  set  of  cracks  of  intensity 
<yk>,  at  locations,  <ek>.  The  fracture  intensity  is  related 
through  the  concept  of  fracture  hinge  to  the  hinge  spring  const¬ 
ant  k.  The  softening  of  the  structure  is  related  to  the  fracture 
hinge  spring  constant  and  the  flexibility  of  the  structure.  For 
each  crack,  it  is  defined  by  a  sensitivity  number  4.  For  a  beam 
of  rectangular  cross-section 


U  -  3*Cl-v2>|[  Uf(X)32dX 


(4.  1 ) 


where  El  is  the  flexural  stiffness,  b  is  the  half-depth,  y*a/b 
defines  the  crack  depth  and  f (X>is  the  dimensionless  fracture 
stress  intensity  factor.  The  direct  and  the  inverse  problems, 
relating  the  relative  modal  frequency  changes  <R^>  to  the  frac¬ 
ture  damage  characteristics  < efe ,  4^3  are  thus: 


R.  «=  f .  <e.  ,4.  )  (4.2) 

i  l  k*  k 

(e,4)k  -  g|t<Ri>  <4.3) 

where  R  "l-di./tt.  ),  *>,  and  are  the  pre-  and  post-damage  fre¬ 
quencies  of  the  i-th  mode.  For  k -number  of  cracks,  it  is  neces¬ 
sary  to  have  2k  number  of  measures  of  relative  frequency  changes 
<Ri ) .  Because  of  the  non-linear  nature  of  the  damage  functions. 


t 
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usually  at  least  on*  additional  measure  of  frequency  change  is 
needed  to  detereine  one  unique  solution  fro*  the  set  of  Multiple 
solutions. 

r 

4.2  Experiment  on  the  Validity  of  "Fracture  Hi nneV Modeling 
of  a  Cracked  Section: 

The  experiment  used  twenty  samples  of  aluminum  cantilever 
beam  with  milled  slits  to  simulate  cracks  of  controllable  confi¬ 
guration.  There  are  five  different  single— crack  locations  along 
the  beam  in  the  length  modulated  dimensionless  locations  <e>  « 
<0.01388,  0.250,  0.51368,  0.70833,  0.791661.  Two  crack  depths 
<10,  20)  mils  correspond  to  <y>  « <0.21333,  0.42666),  or  <■*>  * 
<0.006,  0.024).  The  elastic  effect  at  the  cantilever  built-in 
end,  however  small,  is  taken  into  consideration  in  the  analyti¬ 
cal  model  with  an  elastic  torsional  spring.  The  spring  constant 
t  is  determined  experimentally  through  comparison  of  the  experimen¬ 

tal  frequencies  and  theoretical  ones  for  the  undamaged  beam.  The 
elastic  built-in  end,  therefore,  can  be  assigned  a  sensitivity 

number  4  .  The  cracked  beam  is  therefore  modeled  by  a  double- 
o 

hinge  beam.  Using  Equation  <4. 2),  the  relative  frequency  chan¬ 
ges  can  be  plotted  analytically  as  functions  of  the  crack  loca¬ 
tion  for  the  given  crack  geometry.  Figures  4.1  and  4.2  show  res¬ 
pectively  the  Mode  1  and  the  Mode  2  frequency  changes  comparing 
the  experimental  measurements  of  relative  frequency  changes  with 
the  theoretical  curves  for  two  given  crack  geometries.  In  the  fi¬ 
gures,  the  solid  lines  represents  the  theoretical  relative  modal 
frequency  variation,  based  on  the  average  sensitivity  number  at 
the  built-in  end.  The  dash  lines  and  data  points  denote  the  act¬ 
ual  measured  changes  using  individual  built-in  end  sensitivity 
numbers  of  the  samples.  The  accuracy  is  better  than  0.3%  for  the 
beams  with  the  slit  closer  to  the  built-in  end.  For  beams  with 
»  the  slit  closer  to  the  free  end,  the  errors  is  as  high  as  0.7%. 
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Fig. 4.1  Modified  Mode  1  Fre-  Fig. 4. 2  Modified  Mode  2  Fre¬ 
quency  Variations  quency  Variation 

4.3  Experimental Diagnosis  of  Fracture  Oamaoet 

When  the  built-in  end  sensitivity  numbers  of  the  Individual 
sample  beams  are  used  in  the  double  hinge  model  for  the  diagno¬ 
sis  of  fracture  damage,  using  the  measured  modal  frequencies, 
the  locations  <e>  of  the  damage  intensities,  in  terms  of  indivi¬ 
dual  fracture  damage  sensitivity  numbers,  are  computed  by  the 
computer  program  F1S0L,  developed  for  the  diagnostic  problem 
C301,  as  listed  in  Table  4.1.  The  diagnostic  solution  can  also 
be  expressed  in  the  <e,4>  phase  plane  as  shown  typically  in 
Figures  4.3  through  4.7,  representing  five  different  fracture 
positions.  In  ass  five  figure,  solid  circles  denote  curves  for 
the  relative  frequency  changes  of  the  first  mode;  soled  diamonds 
those  of  the  second  mode  and  solid  triangles  those  of  the  third 
mode.  The  intersections  of  all  three  curves  identify  uniquely 
the  solutions  of  the  damage  characteristics  pairs  <e,4>.  The 
diagnosed  damage  characteristics  for  the  samples  were  compared 
with  the  actual  damage  character istics  of  individual  samples. 

The  errors  in  location  close  to  the  built-in  end  are  less  than 
0.05%.  Up  to  around  1%  near  the  free  end. 
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Sensitivity  timber  (B)  («  14  E-3) 


TABLE  4.1 
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Fig. 4. 4  o-4  Phase  Plane 
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The  rtMirch  uitt  thf  experimental ly  obtained  pro-  and  post 
Matures  modal  frequencies  to  provide  a  positive  identification 
of  the  fracture  damage  CIS, 19,22,26,323.  The  diagnostic  accur¬ 
acy  for  the  two-  hinge  model  (one  hinge  to  simulate  the  elastic 
built-in  end  and  another  to  simulate  the  crack)  is  better  than 
1.02  in  crack  locations  and  around  42  in  damage  intensities. 

For  cracks  at  sections  of  low  resisting  moments,  the  frequency 
measures  will  be  at  the  level  of  noise,  resulting  in  large  eri — 
ors.  Table  4.1  shows  poor  measurements  of  frequency  variations 
for  Samples  Nos.  5,  10  and  14.  These  samples  do  have  slits  near 
the  free  end,  where  the  moment  is  small.  There  are,  however, 
other  insufficient  measurements ,  for  which  the  frequency  varia¬ 
tions  are  in  the  measurement  error  range.  As  in  the  design  of 
the  specimens,  the  slits  locations  were  designed  to  avoid  the 
theoretical  inflection  points.  However,  because  of  the  compli¬ 
ance  of  the  built-in  ends,  slits  at  0. 250  and  0.51389  are  in  the 
neighborhood  of  mode-2  and  -3  inflection  points  respectively. 
Yet,  if  the  frequency  variations  in  other  modes  are  measurable, 
insignif icant  change  in  a  particular  modal  frequency  is  an  indi¬ 
cator  that  the  crack  may  occur  near  that  modal  inflection  point; 
while  the  other  modal  frequency  variations  will  serve  to  detei — 
mine  the  damage  intensity. 
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5.0  UNCERTAINTY  THEORY! 

In  the  expression  of  the  d«Mgt  function,  equation  (4.3), 
it  is  presuned  that,  for  k -number  of  cracks,  (2k+l>  erasures  of 
frequency  changes  are  available  for  diagnosing  the  cracks. 
However,  the  n ueber  of  cracks  is  generally  not  known  a  priori. 
Furthereore,  there  eay  not  be  (2k+l)  Measures  of  frequency 
changes  available.  Uncertainties  thus  result  associated 
essentially  with  a  daeage  configuration  with  Multiple  cracks. 
The  lack  of  a  priori  knowledge  in  excitation  could  also  lead  to 
ei si nf creation  in  Modal  damage,  leading  toward  large  errors  in 
diagnosis.  The  investigation  addresses  these  uncertainties. 


5.1  Frequency  Cross-Overt 

The  uncertainty  of  frequency  cross-over  can  be  illustrated 
with  two  Modal  frequencies,  and  (without  loss  of  genrality 
let  at  undaeaged  state.  After  cracking  occurs,  the  new 

Modal  frequencies  are  correspondingly  and  Zy  If  the  crack 
occurs  near  the  inflection  point  of  the  i-th  nodal  shape,  but  if 
the  location  happens  to  be  near  the  Max i mum  nonent  section  of 
the  j-th  nodal  shape,  nay  not  differ  too  Much  fora  w^.  It  is 
conceivable  that,  for  the  new  nodal  frequencies,  we  nay  record 

In  that  case,  without  the  knowledge  of  the  actual  damage 
conf iguration,  the  diagnostic  assenblage  of  frequencies  nay  well 
confuse  to  be  the  new  frequency  of  the  i-th  node  and  to  be 
the  new  frequency  of  the  j-th  node.  The  phenonenon  is  a  frequen¬ 
cy  cross— over,  which  would  commonly  occur  for  adjacent  Modal  fre¬ 
quencies  of  close  Magnitudes.  The  diagnostic  result  nay  be  dis¬ 
astrous.  However,  with  the  use  of  transmissibility  theory,  for 
which  the  structure  can  be  excited  at  chosen  frequencies,  the 
problem  of  frequency  cross-over  can  be  alleviated. 
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5.2  Close-Facked  Multiple  Crackst 

The  investigation  indicated  that  the  first  few  characteris¬ 
tic  values  of  a  beaa  with  closely  spaced  Multiple  cracks  are*  in 
general  close  to  those  of  a  beae  with  an  equivalent  single  crack 
whose  sensitivity  number  is  approx isatel y  equal  to  the  sue  of 
the  individual  sensitivity  nusbers  of  the  cracks  on  the  original 
beaa.  Furthermore,  the  location  of  the  equivalent  crack  is  gen¬ 
erally  within  the  region  where  the  group  of  cracks  is  located. 
Hence,  equivalence  of  closely  spaced  cracks  to  a  single  crack 
iaplies  that,  in  the  process  of  daaage  diagnosis,  it  is  iapossi- 
ble  to  distinguish  between  closely  spaced  Multiple  cracks  and  a 
single  crack. C 19,26,303  The  lieit  spacing  for  a  cantilever  beam 
is  shown  in  Fig.  5.1. 


Fig.  5.1  Limit  of  Spacing  for  Multiple  Cracks 
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5.3  ifacertainty.  from  Inadequate  Heasureme nti 

The  relative  frequency  changes  at  1  otter  modes,  for  a  struc¬ 
ture  with  eultiple  cracks,  ey  be  the  ease  or  within  the  measure¬ 
ment  tolerance  of  a  structure  with  fewer  cracks.  The  clarifica¬ 
tion  can  only  be  brought  about  with  sore  Measurements  of  fre¬ 
quency  changes.  Akgun  and  Ju  £26,301  illustrated  with  the  com¬ 
parison  of  a  single  major  crack  against  a  group  of  ten  uniformly 
spaced  minor  cracks.  The  first  two  modes  yield  the  same  frequen¬ 
cy  change.  Even  the  third  mode  shows  only  a  difference  of  10~*v 
which  is,  for  all  practical  purposes,  unmeasurable. 


For  illustration,  the  following  relative  frequency  changes 
are  assumed  to  have  been  computed  from  the  measured  frequencies 
of  a  damaged  cantilever  beami  Rj*. 0579,  R^-,0586,  ^*.0591.  It 
is  to  be  determined  whether  the  damage  mainly  consists  of  one 
crack.  The  (actual)  characteristic  values  of  a  damaged  beam  are 
computed  from 


<j> 


j(J>  - 

ac  pu  y*^j* 


(5.1) 


where  ($u'"  >  are  the  character istic  values  of  the  corresponding 
undamaged  beam,  the  first  three  of  which  for  a  cantilever  are 
<1 .6751 ,  4.6941,  7.8348>.  From  Equation  (5.1),  the  third  modal 

ged  cantilever  beam  would  be 


characteristic  value  of  the  d 

j!?,«7.6192.  The  first  two  actual  characteristic  values  yield 


ac 


(3). 


one  solution  for  a  single  crack  located  at  e^  *.372  with 
0.128  (corresponding  to  relative  crack  depth  of  y>.5S  for  a  slen¬ 
derness  ratio  of  b/L  ■  .05).  These  equivalent  values  are  then 
used  in  Equation  (4.2)  with  k«l  to  determine  the  third  character 
istic  value  due  to  the  equivalent  crack.  The  computation  thus 

Then,  the  error  is 


yield  ri^-7.6176. 

•9 


error  -  (*<?>  -  *<3>)/P<3> 


.0002. 


ac  req  ac 
It  is  quite  reasonable  to  accept  a  solution,  whose  erri 


is  less 
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than  O.iX.  However ,  the  given  R^  data  were  actually  ganaratad 
*ro<*  a  problee  with  tan  cracka  o f  equal  intensity  CP  was  taken 
to  be  .01  (i.e.  *  y».2)  -for  all  the  cracks!  located  such  that 

e^.094  for  i  ■  2,3,..  ,10,  where  ei  defines  the  spacing 
between  cracks.  Thus,  in  this  case,  even  though  the  actual  and 
equivalent  values  satch  closely,  a  firs  diagnosis  cannot  be 

reached. 


6.0  PROBABILISTIC  DISTRIBUTION  OF  MULTIPLE  CRACKS  IN 


STRUCTURES* 

The  investigation  considers  the  development  of  the  probabi¬ 
listic  methodology  for  the  prediction  of  multiple  cracks  distri¬ 
bution  in  a  structure  of  beam  elements  associated  with  indivi¬ 
dual  modal  oscillations.  The  probabilistic  Mature  of  crack  dis 
tribution  can  then  be  used  for  the  probabi 1 istic  diagnosis  of 
crack  damage  (depth)  and  its  location  (spacing)  under  random 
loading  and  resolve  some  of  the  intrinsic  uncertainties  in  the 
modal  theories  of  fracture  diagnosis.  The  structural  system  con 
siders  some  randomness  of  material  strength.  The  arresting  frac 
tore  toughness  is  characterized  as  a  random  variable  with  the 
appropriate  probability  distribution.  The  application  of  LEFM 
in  connection  with  the  stress  relief  effect  due  to  the  presence 
of  crack  suggests  a  means  of  predicting  depth  and  spacing  of  ten 
sion  cracks  at  a  given  random  modal  oscillation.  The  resulting 
redistributed  random  bending  stress  (moment)  distribution  is  a 
measure  to  compute  the  subsequent  crack  state.  With  postulation 
that  secondary  cracking  is  dominantly  affected  by  its  immediate¬ 
ly  preceding  crack,  the  process  of  the  successive  cracking  can 
be  treated  as  a  Harkov  process.  The  analyses  are  performed, 
under  these  probabilistic  assumptions,  for  the  first  few  rep¬ 
resentative  normal  modes  of  interest.  The  probability  distri¬ 
bution  of  the  overall  structural  system  therefore  is  obtained 
depending  ext  weight  distribution  of  modes  for  a  particular 
excitation  spectrum  C283. 


6.1  Crack  Deptht 

To  estimate  the  depth  of  tension  edge  cracks  in  a  beam,  a 
fundamental  fracture  mechanics  theory  is  required  to  assess  the 
entire  problem  involved  in  crack  initiation,  propagation,  and 


-22- 


arresting.  For  simplicity  in  the  present  study,  only  model  I 
crack  is  considered  to  initiate  if  the  bending  stress  at  the 
extreme  surface  exceeds  the  cracking  strength,  or.  The  crack 
propagates  when  the  crack  edge  intensity  factor  Kj  reaches  a 
critical  value  Kjc,  fracture  toughness.  For  crack  arresting, 
we  will  assume  that  a  crack  will  continue  to  propagate  as  long 
as  Kj  remains  greater  than  some  critical  value  Kjo  denoted  as 
the  arresting  fracture  toughness.  Kjc  My  be  obtained  via  well- 
defined  test  procedures.  However,  there  is  no  test  data  avail¬ 
able  for  KIo>  Since  both  values  of  KJc  and  Kjq  are  expected  to 
be  of  the  same  order  of  Mgnitude,  it  does  not  seem  unreasonable 
to  treat  Kjd  as  a  uniform  random  variable  in  (0,Kjc>. 

So  long  as  the  rate  of  crack  growth  does  not  approach  the 
speed  of  propagation  of  an  elastic  disturbance,  the  computed  Kj 
on  a  basis  of  static  elastic  theory  should  be  a  good  approxima¬ 
tion  to  describe  the  arresting  process.  Analytical  solutions  for 
Kj  only  exist  for  selected  simple  cases.  The  solution  of  an  edge 
crack  in  half-plane  is  approxiMtely  applied  to  our  beam  sx>del 
with  edge  crack  length  d  in  the  presence  of  the  bending  stress 
field. 

K.  «  f  1 .  lZ/xd  - * o.  <l-2y/t>dy  <6.1> 

1  jO  *<d  -y*> 

where  is  the  outer  fibre  bending  stress.  Kj  is  readily  ob¬ 
tained  by  integrations 

c <k_/o.  )  «  <*/4>n  -  n3  -  f<n>,  <6. 2> 

X  b 

where  c  *  /*7t/4.48,  n  ■  d/t,  and  t  is  the  depth  of  the  beam. 


The  variation  of  Kj  with  crack  depth  parameter  f)  is  schema¬ 
tically  plotted  in  Figure  6.1,  wherein  f  (t»  assumes  its  Mximun 

value  f  at  tj  .As  the  figures  shows,  if  c  <K.  /a.  >  >f  ,  no 

max  max  Zb  max  ’ 

crack  will  occur.  Thus,  the  cracking  strength  ar  can  be  related 
to  Kjc  by  applying  <6. 2)  at 


C<KIc/ar>  "  'max’ 


(6.3) 
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Fig.  6. i  Variation  of  Kj  with  Crack  Length 

which  yields  ar  ■  cKj  /f^^.  Equivalently,  it  sieply  says  that 

crack  occurs  if  ofa>or.  Also  shown  in  the  saee  Figure,  if  cKi0,ab 

is  less  than  F^gv ,  there  are  two  positive  real  value  solutions 

r)  «nd  T)  with  q  <q  .  q  is  interpreted  as  the  normalized  cri- 
c  o  c  o  c 

tical  crack  length  For  unstable  cracking  to  propagate  and  qQ  way 
be  treated  as  the  Final  crack  length  being  arrested.  Since 
is  assumed  to  vary  within  (0,Kjc>,  the  prediction  oF  q  needs  the 
random  variable  analysis.  The  normalized  mean  crack  depth  Etd/tl 
can  be  computed  asi 


ECd/tl  - 


r O. 57 (z/3) <0./0  > ,  0K<0„ 

'  b  r  b  r  C6.4) 

l  CO. 52-0. 15c os C54/3) -0.5cosF-0. 75c os (4/3) 3 (z/3) co*4 


ab>ar 


where  cos-*  «  -e  /a.. 

r  D 


6.2  Crack  Spacing  and  Stress  Relief* 


The  quantity  oF  principal  importance  to  the  problem  of 
crack  spacing  is  the  stress  perturbation  in  the  plane  in  which 
cracks  originate,  usually  the  Free  surface.  The  elastic  stress 
perturbation  due  to  the  presence  of  crack  in  the  semi-infinite 
medium,  used  with  the  modified  Griffith  theory  of  macroscopic 
Fracture,  suggests  a  means  of  predicting  spacing  of  tension 
cracks  in  terms  of  stress  Field  and  measurable  properties.  Since 
the  modal  theories  of  Fracture  diagnosis  are  non-destructive  and 
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ts  are  carried  out  in 


us*  Ion  level  oscillation,  the  eeasureaen 
the  elastic  range  of  the  structure,  within  which  the  characteri¬ 
sation  of  the  fracture  daaage  can  be  defined  with  LEFH.  Figure 
6.2  illustrates  the  noraal ized  surface  tangential  stress  relief, 
ge<?)  and  gfi(f),  respectively,  for  step  and  linear  distribution 
of  noraal  stresses  on  the  crack  surface.  The  non-dieensianal 
variable  ?  denotes  the  distance  aeasure  froa  the  crack. 


Fig. 6. 2  Noraal ized  Edge  Surface  Bending  Stress  Relief 

For  a  beam  with  crack  depth  d  and  the  linear  initial  bend¬ 
ing  stress  distribution  across  the  thickness  t,  the  bending 
stress  relief  g(f>  can  be  computed  froa  the  combination  of  gft(?> 
and  g^if)  within  the  linear  elastic  range. 

g (f )  «  g  <f>  -  2(d/t>  g_<f>.  <6.5> 

®  e 

6.3  The  Nodal  Cracking  Probability  of  Beam  with  Uniformly- 
Distributed  Cracking  Strength  Subjected  to  Random 
Amplitude  Sinusoidal  Loading* 

One  commonly  used  method  of  numerical  simulation  of  earth¬ 
quakes  is  to  expand  a  periodic  function  into  a  series  of  sinu¬ 
soidal  waves. 
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(6.6) 


n 

6 (t )  *  £  A.  sin(w.t+4), 

i-1  1  1  x 

where  A^  is  the  amplitude,  is  the  excitation  frequency,  and 
4^  is  the  phase  angle,  respectively,  of  the  i-th  contributing 
sinusoid.  Introducing  randomness  to  4^ ,  a  wide  band  random  ex¬ 
citation  can  be  generated.  In  the  present  paper,  a  simple  sinu¬ 
soidal  ground  motion  6(t>  *  A  sin(ut)  is  taken  with  A  as  random 
variable. 

For  clarity  convenience,  the  notation  *(j,k)  beam"  is  used 
to  identi-fy  the  state  of  a  beam  being  in  j-th  mode  and  possess¬ 
ing  k  numbers  of  cracks.  The  associated  quantities  are  then  exp¬ 
ressed  by  subscript  state  parameter  (j,k).  Thus,  for  a  general 
N-crack  beam  system  (d^k>  e^k>  is  the  fracture  damage  charactei — 
istic  pair  denoting  the  k-th  crack's  depth  and  its  location  at 
j-th  mode. 

The  spacing  s^k  is  defined  as  the  location  difference  bet¬ 
ween  two  sequential  cracks,  e..  -  e.  .  The  expression  for 

the  random  edge  surface  stress  S^k(f,t),  its  associated  positive 
quantity  X..(?tt)f  and  its  related  positive  random  quantity 

jk 

Yjk<?,t)  beyond  the  random  crack  strength  R  are: 

Sjk<f,t>  «  A  sin(«t)mjk<D  ,  (6.7) 

Xjk(*,t)  «  Sjk(f,t)  HCSjk(f,t)3,  (6.8) 

and 

Yjk(f,t)  -  ISik<e,t)-RJ  HCSjk(e,t)-R3,  (6.9) 

where  A  and  R  are  uniformly  distributed  in,  (0,n or>  and  (0,or> , 
respectively,  H  is  the  Heavyside  function,  and  p^k(C)  i*  the 
corresponding  modal  surface  stress  shape  of  the  (j,k)  beam. 
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6.4 


The  solution  o-f  the  enn  time  cycle  o-f  cracking,  the  proba¬ 
bility  of  cracking,  and  the  axial  probability  distribution  o-f 
cracking  can  be  -forealized  with  procedure  steps.  Steps  1-3. 

Steps  4  and  S  compute  the  mean  spacing  between  sequential  cracks 
and  the  mean  depth. 


Step  <  1 )  For  (j,k)  beam  with  damage  pair  (d^^e^)  ,  update  the 
modal  bending  shape  -function  and  obtain  Sj^ffyt),  Xjk<P,t) 

and  Yjk(?,t)  -from  Equations  (6. 7-6. 9). 

•jk'1' -  •*«<»  <6'i0> 

Step  (2)  Linder  assumption  2,  the  mean  time  cycle  m  -for  the  suc¬ 
ceeding  cracking  and  its  cracking  probability  F.  .  , (£ _ ,mT>  can 

9  K+ 1  MX 

be  evaluated  in  conjunction  with  X^  ,t)  t 


Xjk<fmax»t) 


^jk^max’*”*3 


FrCA  sin(wt)p..(f  ,t)  >R3 

J  K  RloX 


1  -  1/  <n  sin<«t)pit,(f^v)  ,  l<npjkCf(njix)sin(ut> 


jk  max 

otherwise 


,pI-f 


Xjk't«x’T,dT3  ' 


and 


/<!  -  expt 


Fj,k+l(f.ax‘«T) 


-J1 


*jk<fmax’T,dT3>’ 


*  1 


-expt-f 


X..  <f  ,t)dtl. 
,  j  k  max  * 


(6.11) 


<6. 12) 


(6.13) 


Step  (3)  With  the  evaluation  o-f  <  EtY^djt)!  >,  the  axial  pro¬ 
bability  density  -function  p.  .  ,(f,t)  can  be  computed! 

J  ,  K  +  l 

EtVjk<t,t>3  -  Cn  .in<.t),jk(f)/6-  \  *Jk«t>3, 

(6.  14) 
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<ECY.k<f,t)3>-=i  JoEtVjlt<t,T»  dr 

-  2  LEtYjk(t'T'  dT‘ 


(6.  IS) 


•nd 


<ECY  <f,«T>l> 

■f-jfcCV-*k<e*“T>3> 

<EtY.k(f  ,T)3> 

Ji^CYik‘t'T,I>  dt 


(6k  16) 


where  e^k<£<l. 


Step  (4)  Computations  -for  mean  spacing  s.  ....  and  crack  location 

i  •  k+l 


j  ,k+l 


are: 


*J,k-l  "  Jej*"*jk>pj,k+l<f »*T>df- 


(6.17) 


and  e.  ...  «  s .  ,  . ,  *  e... 
J « k+1  j,k+l  jk 


Step  (5)  For  the  (k+i)-th  crack  produced  at  e.  ...»  the  temporal 

J  9 

mean  stress  <ECX  (e  . ,bT)3>  up  to  time  mT  ist 

JK  Jy 


<ECXjk<«J,k*l’"T)3> 


mT 


r- 

*  c\ 


*  j. 


“^k^j.k+i*^  ^  * 


c  X,.  (e 


Jk  j,k+l 


|T)  dr. 


(6.18) 


Replacing  eb  in  Equation  (6.4)  by  <ECX^k(e^  k+l»mT)D>  yields 
depth  d^  k+J.  This  new  damage  pair  (d^  k+l»  «j  k+1>  defines  a 
new  crack  state,  namely,  (J,k+1>  beam.  The  loop  completes  and 
returns  to  step  (1). 
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7.0  MODELING  OF  ELASTIC-PLASTIC  8TRUCTllRESt 


This  investigation  established  a  procedure  -for  the  probabi¬ 
listic  characterization  of  a  daaaged  structure*  which  is  modeled 
as  a  nonlinear  elastic-plastic  MDF  system.  The  random  excitation 
may  be  either  stationary  or  nonstationary.  The  stiffness  matrix 
is  nonlinear  to  simulate  the  elastoplastic  behavior  of  a  damaged 
structure.  The  stiffness  matrix  is  also  random  to  characterize 
the  material  and  environmental  variations.  The  governing  stoch¬ 
astic  differential  equation  is  resolved  into  one  for  the  mean 
response  and  another  for  its  random  component.  Responses,  their 
statistical  moments  and  cross-moments  are  solved  with  discrete¬ 
time  recurrence  formulations.  The  probability  of  structural 
damage  or  the  structural  reliability*  is  then  estimated  by  the 
upper  bound  of  the  cumulative  energy  dissipation.  The  formalism 
of  approach  in  formulating  the  solutions  for  a  generic  class  of 
MDF  nonlinear  problems  with  Prandt-Reuss  material  permits  ready 
adaptation  to  REM  analysis.  117,313 

7.1  Formulation; 


The  governing  differential  equation  of  motion  for  a 
non-linear  structural  framework,  modeled  as  a  discrete  MDF 
system*  can  be  written  as 


mz  ♦  cz  ♦  k(z)  z  *  f 


(7.1) 


where  m*  c*  k<z>  are  the  NxN  mass,  damping,  and  stiffness  mat¬ 
rices,  f  is  the  external  load  vector \  z«z,z  are  the  displacement, 
velocity,  and  acceleration  vectors  of  the  system.  In  this  inves¬ 
tigation  m,c  are  assumed  to  be  deterministic  and  constant.  The 
stiffness  is  represented  by  a  matrix  of  random  variable  which  may 
correlate  to  the  response  z<t>.  The  quantities  f*  lc(x>,  and  z 
can  be  resolved  as  follows: 
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f  «  *+F,  k(z>  -  X(p)+K«  z  «  C+Z ,  (7.2) 

•there  9 ,  X  end  P  ere  Men  velues  thet 

ECf3  «  9,  ECk (z> 3  -  X(p>,  ECz3  -  p.  (7.3) 

It  is  noted  that,  (7.2),  the  nonlineer  properties  of  stiffness 
ere  essumed  to  be  reflected  by  its  mean  component.  The  random- 
ness  of  stiffness  is  represented  by  K  which  is  e  eetrix  of  rendom 
veriebles.  F  is  e  non-stationery  bend-limited  white  noise.  Also 
it  is  noted  thet  the  rendom  quantities,  F ,  K  end  Z  introduced  in 
(7.2)  ere  ell  zero-mean.  Substitution  of  (7.2)  into  (7.1)  yields 

m  »«  mm 

m(p+Z)  +  c(p+Z>  +  EX(p)+K3 (p+Z>  *  *+F.  (7.4) 

Taking  expectation  on  both  sides  of  the  above  equation  results  in 
the  mean  values  equation. 


•  •  e 

mp  +  cp  +  X(p)p  *=  p-ECkz  3 .  (7.5) 

The  difference  of  (7.5)  end  (7.4)  yields  the  equation  for  the 
random  component. 


•  t  • 

mZ  +  cZ  +  X(p)Z  *  F  +  ECKZ3  -  KZ  -  Kp.  (7.6) 

In  (7.5)  end  (7.6),  KZ  end  ECKZ3  being  the  product  of  two  rendom 
quantities,  ere  higher  order  terms.  By  neglecting  these  terms, 
(7.5)  end  (7.6)  ere  reduced  to  the  following  expressions! 


m  m 

mp  +  cp  *  R<p)  »  9,  (7.7) 


mZ  ♦  cZ  ♦  X(p)Z  -  F  -  Kp.  (7.8) 
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The  omission  o-f  the  higher  order  terms  ere  postulated  -for  the  HDF 
systems.  For  SDF  systems,  the  errors  in  response  and  its  statis¬ 
tics  resulting  from  omission  of  the  higher  order  terms  is  refer¬ 
red  to  £17,313. 


7.2  Nonlinear  Models 

The  investigation  treated  a  generic  class  of  elasto-plastic 
material  that  satisfies  the  Prandtl-Reuss  relationship.  For  a 
beam  with  symmetric  cross  section,  the  moment  for  any  section  at 
location  x,  measured  longitudinally  along  the  beam,  at  pure 
'bending,  is 

M(x>  *  Ely"  -  E  f  c  <x,h)  hw(h>  dh,  <7. 9) 

where  E  is  the  Young's  modulus,  I  is  the  moment  of  inertia  of  the 
cross  section,  y"  is  the  second  derivative  of  the  deflection  with 
respect  to  the  coordinate  x,  c  is  the  half  depth  of  the  symmetric 
cross  section,  eQ<x,h>  is  the  permanent  set  at  locations  x  and 
h.  All  external  load,  without  loss  of  generality,  are  resolved 
at  nodal  point.  Hence 


X  c 

Ely'  -  E  J  du  |  «Q<u,h)hw(h>dh  +  ^CjX2  ♦  C2x  ♦  C3,  C7.10) 
X  V  c 

Ely  «  eJ  dv  J  du  |  *o<u,h>hw(h>dh  +  Jc^x3*  2C2*2*  C3X  +  C4 

OOC  .  -m  ... 


(7.11) 


The  constants  ,  i  ■  l,..., 4  are  prescribed  from  boundary  condi¬ 
tions.  The  axial  strain  can  be  evaluated  with  an  analogous  ap¬ 
proach. 


<x> 


EA 


[P  +  E  «  <w,h>  w (h)  dh], 


C7.12) 
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where  P  is  the  axial  force,  L  is  the  length  of  the  beam,  A  is  the 
cross-sectional  area.  The  permanent  sets  are  yet  to  be  deter — 
mined.  For  that  purpose  an  iteration  scheme  is  developed  using  a 
finite  difference  method  that  the  displacement  at  time  t^+j  can 
be  evaluated  from  (7.7)  by  using  the  central  difference  method; 

*Vl  “  A~1C2mpj  +  +  (7. 13) 

where,  A1  *  ~  cflt  ♦  m,  and  A3  *  2  cAt  ~  •“  (7.14) 

It  is  assumed  that  the  system  starts  at  rest. 


In  the  iteration  scheme,  the  neutral  axis  does  not  change 
throughout  the  computation.  Because  of  this  assumption,  the 
permanent  set  will  first  converge,  then  alternate  between  two 
values.  In  such  cases  an  approximation  for  the  permanent  sets 
can  be  established  by  averaging  the  two  values. 


Since  (7.6)  characterizes  the  random  component  of  the  struc 
tural  response,  the  displacement  response  at  time  tJ+l  can  be 
solved  by  using  central  difference  approximation.  Namely, 

zj.i  *  fl7‘  ‘«2jzj  *  Vj-i  ♦  ct2fj  ~  at2KV  ■  ,7- l5) 

where  Aj  and  A-j  are  given  by  (7.14)  and 

A2J  «  2m  -  (7.16) 


in  which  X(p^)  is  the  equivalent  stiffness  matrix  at  time  t^;  the 
s-th  column  in  X(/ij>  is  by  definition. 


8R(p.) 


(7.17) 


The  response  covariance  matrix  can  then  be 

established. 
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7.2 


Damage  Diagnosis* 

Structural  damage  resulting  fro*  lout-cycle  loading  may  use 
the  cumulative  cycle  peraanent  met  of  the  cyclic  energy 
dissipation,  if  it  is  postulated  that  the  cumulative  damage  is 
related  to  the  energy  dissipation  in  the  system.  In  a  small 
element  of  a  beam,  the  energy  dissipated  in  the  system  due  to 
material  nonlinearity  can  be  expressed  as 


E  *  flV  E  f€C*(t>  -  c  (t  >  3  da.  (7.18) 

s  jo  o 

The  total  energy  dissipation  for  a  beam  due  to  material  nonlinear — 
ity  is  therefore  obtained  by  summing  up  all  the  small  elements, 
resulting  in 


«  E  4V  E  Etc  (jilt)— *(  jilt)  3  tc  « j+l>At)-* ( jAt)  3  (7.19) 

ilV  ° 

where  Ey  is  the  total  energy  dissipated  in  a  beam.  It  is  noted 
that  e(t),  eo(t)  *re  obtained  from  (7.7)  which  represents  the 
mean  of  the  response.  In  this  sense,  the  energy  dissipation  cal¬ 
culated  by  (7.19)  can  be  used  to  represent  the  mean  value  of  the 
energy  dissipation.  The  result  is  used  to  predict  an  upper  bound 
of  the  energy  dissipation  areasure,  considering  Markov  inequality, 

Pr fx >d>  1  (7.20) 

where  d  is  a  positive  constant.  If  total  energy  dissipation  is  a 
maximum  damage  criterion,  the  probability  of  damage  is  therefore 
defined  by  (7.20)  with  d  representing  the  material  constant  of 
permissible  total  energy  dissipation.  Numerical  examples  are 
referred  to  £17,313. 
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8.0  RELIABILITY  OF  STRUCTURES  WITH  STIFFNESS  AND  8TRENBTH 

DESRADATIONt 

When  a  structure  of  frictional  materials,  such  as  reinforced 
concrete,  is  subjected  to  strong  randoe  excitations,  the  struc¬ 
ture  eay  undergo  inelastic  deformations  during  certain  cycles  of 
loading,  with  associated  cyclic  degradation  in  stiffness  or 
strength,  or  both.  The  exact  nature  of  system  degradation  is  a 
function  of  the  structural  materials  and  the  conf iguration ,  and 
may  vary  considerably  from  structure  to  structure.  Basically, 
the  deteriorated  phenomenon  is  due  to  the  extension  of  crack  in 
the  concrete,  bond  deterioration,  bar  slippage,  shear  deforma¬ 
tion,  and  inelastic  deformation  of  reinforcement.  The  most  im¬ 
portant  factor  is  the  opening  and  closing  of  the  crack  in  the 
structure  that  alternate  between  compression  and  tension  during 
the  response  cycle.  Thus,  the  opening  and  closing  of  these 
cracks  may  eventually  lead  to  a  deteriorated  stiffness  and 
strength.  As  a  consequence,  energy  is  dissipated  through  the 
degradated  hysteresis.  It  has  been  experimentally  verified  by  Ju 
et  al  E14D  that  the  rate  of  degradation  is  related  to  the  energy 
dissipation  through  the  degrading  restoring  hysteretic  loop. 

Quite  frequently,  the  excitations  of  the  structural  system 
are  not  predictable.  Examples  of  such  loading  sources  are  earth¬ 
quake,  wind,  aerodynamic  loads,  etc.  Under  these  types  of  excita¬ 
tions,  the  structural  response  apparently  behaves  randomly.  The 
random  characteristics  of  the  loading  sources  together  with  the 
random  system  degradation  lead  to  the  desire  of  a  coherent  damage 
model  that  can  be  used  to  assess  the  reliability  of  such  a  sys¬ 
tem.  In  practice,  the  randomness  of  the  excitation,  together 
with  the  randomness  of  the  system  degradation,  will  lead  to  the 
randomness  of  the  deteriorated  restoring  force  to  be  a  random 
process,  which  we  refer  as  structural  random  noise.  The  predic¬ 
tion  of  reliability  and  assessment  of  damage  depend  upon  the 
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proper  Modeling  of  such  structures,  taking  into  consideration  the 
random  characteristics  of  the  Materials  as  well  as  the  excita¬ 
tions.  Therefore,  the  present  investigation  establishes  such  a 
Model  that  can  be  used  to  predict  the  liability  of  a  generic 
nonlinear  structural  systea,  especially  for  those  that  show 
stiffness  and  strength  degradation. 


6. 1  For mu I  at i  on: 

The  nonlinear  systee  to  be  considered  herein  is  a  single- 
degree— of  — freedom  (SDF)  system,  with  the  governing  differential 
equation  of  motion 


•  •  e 

m  z  +  c  z  +  R(z )  «  f<tl,  <8. 1 ) 

where  m,c  are  the  mass  and  damping,  respectively;  fit)  is  the  ex¬ 
ternal  random  excitation,  z,  z,  z  are  the  displacement,  velocity 
and  acceleration  responses  of  the  structural  systea,  respective¬ 
ly.  Rtz)  is  the  system  hysteretic  restoring  force.  In  this  inves 
tigation,  m,  c  are  assumed  to  be  deterministic  and  constant.  Fur 
ther,  it  is  assumed  that  fit)  is  a  band  limited  zero-siean  station 
ary  white  noise  with  constant  power  spectral  density  namely, 

ECf it)fis)3-  pff*<t-s).  <8. 2) 

The  behavior  of  the  system  hysteretic  restoring  force  for  a  sys¬ 
tem  that  shows  stiffness  and  strength  degradation  has  been  stud¬ 
ied  extensively.  Among  those  system-degradation  hysteretic  mo¬ 
dels,  the  Q-hysteresis  can  reproduce  the  behavior  of  system  deg¬ 
radation  in  a  simple  and  efficient  way.  Hence,  it  is  adopted 
here  in  the  present  study  to  describe  the  deterministic  behavior 
of  system  degradation.  The  rules  of  Q-hysteresis  are  summarized 
in  Figure  8.1,  in  which  kj  is  the  initial  stiffness.  The  unload¬ 
ing  stiffness  k^  is  determined  bys 

k_  ■  k,  Cz  /z  3**,  y  ■  0.4.  <8. 3) 

3  1  y  max  * 
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It  is  noted  that  the  parameters  k^,  i«l,...5  in  Figure  8.1  define 
the  rules  of  loading,  unloading  and  loading  reversal .  On  a  simp¬ 
ler  level,  all  these  parameters  can  be  assumed  that  they  do  not 
manifest  any  randomness.  However,  on  a  more  realistic  level,  all 
these  parameters  cannot  accurately  describe  the  original  behavior 
of  system  degradation.  Moreover,  due  to  the  material  randomness, 
even  two  identical  samples  cannot  yield  the  same  hysteretic 
curve.  Therefore,  the  error  that  between  real  behavior  of  system 
degradation  and  the  model  we  developed,  together  with  the  random 
characteristics  of  the  material,  lead  us  to  the  realization  that 
the  hysteretic  restoring  force  must  behave  randomly.  Further,  in 
view  of  the  fact  that  the  restoring  force  is  a  function  of  the 
random  displacement  response,  the  randomness  of  the  hysteretic 
restoring  force  must  be  time  dependent.  Without  loss  of 
generality,  the  hysteretic  restoring  force  can  be  rewritten  as 

R  ■  Rta^ , . . . ,a^>  (B.4) 

where  ,  i*l,...,5,  c  “z.  By  using  the  Taylor's  expan¬ 

sion,  R  can  be  expanded  about  the  means  of  its  underlying  random 
parameters  to  obtain 
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(8.5) 


R<M. 


•V  +  A  <ai“V  S*  * 


2  E  E  <«*-**  Maj-Mj*  Ja  %a. 

*i Kj  jB^  *  *  J  J  i  j 


R(p 


i*' 


•  «f&>  « 


Mhcrc  ji. 


ECai3,  fi«  i“i , -  - -6.  It  is  noted  that  the  deri¬ 


vatives  in  the  above  equation  are  evaluated  at  ,  i«l,...6.  The 
N(?4t. ..,f.)  is  the  generalized  structural  noise.  The  structural 
noise  arises  -fro®  the  uncertainty  of  the  eater ials,  the  errors  in 
the  aodel  and  the  randoeness  of  the  response.  Apparently «  the 
structural  noise  is  a  wide  band  random  process.  A  typical  such 
structural  noise  is  shown  in  Figure  8.2.  According  to  1173,  the 
structural  noise  can  be  approximately  assumed  to  be  zero-mean 
with  less  than  1  percent  error.  Moreover,  from  Figure  8.2  it  can 
be  seen  that  the  structural  noise  indeed  shows  the  property  of 
zero-mean.  Substitution  of  Equation  (8.5)  into  (6.1)  yields 
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Fig.  8.2  A  Typical  Example  of  Structural  Noise 


(8.6) 


mz  ♦  cz  ♦  R<z>  ♦  N<t)  «  f <t>. 

The  structural  noise  N(t)  is  introduced  in  the  above  equation  and 
is  presented  as  a  function  of  tiee  for  sisplicity.  In  order  to 
specify  the  characteristics  of  a  randoe  process,  the  probab  ility 
structure  of  such  a  random  process  oust  be  given.  Namely,  the  pro¬ 
bability  density  functions  of  the  random  process  up  to  infinite 
order  must  be  known.  In  practice,  this  is  impossible.  However, 
the  difficulty  is  alleviated  by  adding  more  restriction  when  eva¬ 
luating  the  structural  noise.  In  the  present  investigation,  the 
restriction  is  made  by  adding  the  energy  dissipation  as  a  given 
condition.  In  such  a  situation,  the  autocorrelation  function  of 
the  structural  noise  can  be  evaluated.  According  to  the  above 
statements,  if  the  wide  band  random  noise  assumption  is  made,  the 
autocorrelation  function  can  be  written  as 

ECN<t)N(t+s) |Ed,  z(t>3  *  S(t )  ff  <s) . 

It  is  noted  that  S(t)  is  time  dependent  which  reflects  the  nonsta¬ 
tionary  characteristics  of  N<t>. 

Furthermore,  if  the  input  excitation  f (t)  is  assumed  to  a 
band-limited  white  noise,  it  becomes  necessary  that  the  theory  of 
Markov  process  needs  to  be  used  here.  In  view  of  this,  the  tran¬ 
sition  probability  density  function  is  defined  as 

P  «  P<t,C,t+At,y>  -  P<z(t+At>-y|z(t)«C,  Ed>.  (8.8) 

Equation  (6)  implies  that  the  random  process  at  time  t+At  is  eva¬ 
luated  with  the  conditions  that  not  only  z(t)»C  has  to  be  given, 
but  also  the  energy  dissipation  sequence  Ed  must  be  realized. 

The  energy  dissipation  sequence  Ed«  therefore,  may  be  viewed  as 
the  random  environment.  6iven  different  environments,  the  transi¬ 
tion  probability  that  governs  the  random  evolution  varies. Simi lar 
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to  <8.8>,  the  first  end  second  eoeents  of  the  ran don  increments 
that  ere  associated  with  z<t)  can  be  defined  as 


!*“>  ■  aiiS  at  Edi, 

°Z<t’  *  aliS  at  Etto2|  Ed3‘ 


<8.9) 

(8.10) 


The  Kolmogorov  backward  equation  still  holds  if  the  deriva¬ 
tions  are  followed.  The  reason  is  that  the  given  random  environ¬ 
ments  didn't  change  the  characteristics  of  )i<t)t  cr(t)f  and  P. 
Hence,  from  <8.8  -  8.10)(  the  Kolmogorov  backward  equation  is 

_  agi.tiCxt „  ..  apit,c.t;ay>  +  2?  a?P<t1c,*;,y>  (e  ... 

at*  p'*'  ac  2  ac*  <b.id 

where  t'  *  t+At,  with  initial  condition 


P  <t , C ,t '  ,y)  «  tf <C-y) . 


(6.12) 


It  is  reasonable  to  assume  that  P  is  stationary  within  time 
Ct,t'3  if  At  is  small.  Based  upon  this  assumption,  (8.11)  becomes 


.BP<t,.C,t'.,yl  „  +  s£lll  ' ,Y>  13) 

at  ac  2  ac*  * 

2 

The  evaluation  of  p<t)  and  a  (t)  are  discussed  briefly  as 

follows.  Let  Az*=Au+AZ  where  u*ELz  3 .  When  taking  conditional 

expectation  on  Az,  from  Equation  <8.6)  it  is  easy  to  show  that 

_t 

M<t)  =  u <t )  +  JecJ  exp<-jjj<t-T>>f  <T>dr3  *  0,  <8.14) 


since  f<t)  is  zero-mean,  f <t)  and  Erf  are  independent.  Similarly, 
from  Equation  <8. 6) 


a2  (t) 


.t'  .t  ‘ 


Atio  At  ECc_2JtdTJtdV  h(T>h<Y>9<t'-r)V<t'-y) jz<t)«C,Ed>3 


<8.15) 
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in  which  f(t)^lt)-W(t),  h  (t)«l-exp  (-ct/m)  ,  where  uwt  have  baan 
aade  t ha  facts  that  Au*0 (At)  and  f (t>  is  zero-mean. 

Equation  (8.15)  can  ba  raducad  further  by  noting  the  fact 
that  EtN(t)f (t) 3*0  since  f(t)  is  white  noise  C173.  The  integra¬ 
tion  can  be  carried  exit  to  result  in  C213 

w2..  c“2  Cv4f+S(t)3tl-exp<-ct/e)32.  (8.16) 

2 

It  is  noted  that  a  (t)  is  also  a  random  variable  at  tiee  t  since 

it  contains  S(t).  Also  e2(t>  approaches  to  c  2Cg<f<f+S(t)  3  as  t**eo. 

The  solution  of  Equation  (8.13)  can  be  solved  using  the  Fourier 

2 

transform  if  p(t)  and  a  (t)  are  given  by  <8.14  and  8.15),  respec¬ 
tively.  The  solution  is  6aussian  distribution  with  parameters  t 
and  V;  namely, 

P(t,C,t',y,S(t))  «  PCz  <t  *  >*y |z (t)*C,  E  )3  * 

“V2W  ^<y-C)/«/V>  (8.17) 

where 

V=c”2Ce+f +S  (t )  3 (t+^ffit  1-exp  <-jj£)  3+jj®!  1-exp  (-2^)  3> 

It  is  noted  that  the  transition  probability  given  by  (8.17)  con¬ 
tains  a  random  variable.  This  means  that  the  characteristics  of 
the  transition  probability  is  random  and  its  value  depends  upon 
the  given  condition  of  the  energy  dissipation.  This  satisfies 
the  original  assumption  which  states  that  given  different  environ¬ 
ments,  the  transition  probability  that  governs  the  random  evolu¬ 
tion  varies.  The  transition  probability  that  contains  random 
variables  are  called  Markov  chain  in  random  environment  (MCRE). 


8.2  Discrete  Formulation! 

The  transition  probability  established  above  is  based  on  the 
assumption  that  z(t)  is  continuous.  However,  in  digital  compu¬ 
tation,  a  discretized  form  is  necessary.  Also,  the  nonlinear 
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characteristics  of  R(z)  and  the  nonstationary  property  of  N(t) 
•ake  the  transition  probability  solved  above  only  valid  within 
saall  time  intervals.  Hence,  these  conditions  necessitate  the 
use  of  discretized  fora  of  <6.17).  If  discrete  fora  is  used, 
Equation  <8. 17)  becoaes 

<P(n,C,n+l ,y,S  )3  «  P  « 
n  n 

‘  . ’o>  V  ,B-,8> 

It  is  noted  that  zQ,...tzn_j  are  also  put  into  the  given  condi¬ 
tion  since  the  future  aay  not  be  independent  of  the  past  for  this 
case.  The  n— step  transition  probability  is  given  by  (8.20). 

P"”  -  Ed”  ‘  Po  ‘frpj*  <8-19' 

j«l 

where  P^ ,  j«=l,..,n  are  the  transition  probability  at  tiae  step  i, 
and  PQ  is  the  initial  probability  distribution  of  zq.  The  proof 
of  (8.19)  is  referred  to  E213. 


8.3  ftpplicat.ipn: 

There  are  two  cases  that  can  be  considered  here.  The  first 

case  is  that  each  of  the  S  is  independent.  For  such  a  case,  the 

n 

structural  response  is  then  a  random  process  moving  in  the  avei — 
age  environment.  The  mean  of  the  n-step  transition  probability, 
in  this  case,  is  obtained  by  taking  expectation  on  both  sides  of 
(8.19)  yielding 

n 

ECP  3  «  ECP  3  ECP . 3 ,  (8.20) 

°  J«1  J 

where  ECP^3  is  the  mean  transition  probability  at  tiae  j.  This  is 
a  special  case  of  MORE.  It  other  words,  z(t)  is  a  nonstationary 
Harkov  chain  with  one  step  transition  probability  ECP^l. 
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It  is  very  interesting  to  note  that  ECP(n,C,n+l ,y ,8  ) 3  is 
still  s  valid  representation  of  a  probability  density  function, 
where  P(nfC ,n+l ,y,Sn> 1  is  given  by  (8.18).  The  reason  can  be  seen 
by  noting  that 

Etf  P(n,C,n+l,y,S  )  dyl  *  f  ECP(n,C,n+lty,S  )ldy  *  1. 

Jy  n  jy  n  <8.21) 

By  using  the  definition  of  the  expectation,  (8.21)  can  be  re¬ 
written  as 

Jv  dy  Jy  P«".t."+‘.y.sn'  dSn  • 

«  Js  P(Sn)dSnJy  P(nfC,n+l,y,Sn)  dy  *  1,  (8.22) 

where  P(Sn)  is  the  probability  density  function  of  &n  at  tiee 
step  n.  Equation  (6.22)  can  be  used  as  a  tool  for  the  nunerical 
computation  and  will  be  discussed  in  Section  8.2. 


The  second  case  is  that  Bn,  n  «  1,2,...  are  also  Harkovian. 
In  such  a  case,  the  one  step  transition  probability  which  is 
given  by  (8.9).  It  can  be  rewritten  by  using  the  law  of  total 
probabi 1 ity. 


P(j»C, j+l,y,Sj)  « 

-  E  E  P*(j,C,j+l,y,S.)  P(S.|s  >  P(S  )  « 

S.S  J 

J  o 

*  E  E  ..E  P*(j,C,j+l,y,S.)P(S.|s  )..P(S.|s  >P(S  ),  (8.23) 

S  .S  .  ,S  J  J|  J  ' 

J  J-l  o 

where  in  P*(j,f,j+l,y,SJ  is  a  given,  deterministic  value.  If 
we  let 


<P(S.  js^)}  -  l\,  i«l ,. . ,  j  (8.24) 

and 

<P(s  )J  «  T  ,  (8.25) 

o  o’ 


then 


P 


,  j 

<P(  j,C,  j+l,y,S  .)>  «P  (j,C,j+l,y,S.)TTr4-  <8.26) 

J  J  i«l  1 


The  n-step  transition  probability  then  can  be  evaluated  based 
upon  (8.19).  Numerical  examples  are  referred  to  1211. 
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9.0  CHARACTERIZATION  OF  DAMAGE  USING  THE  CONCEPT  OF  STRUCTURAL 


POWER 

This  investigation  aims  to  establish  a  theory,  using  the 
structural  power  to  characterize  the  nonlinear  structures.  Since 
a  damaged  structure  displays  uniquely  the  nonlinear  charactei — 
istic,  the  established  theory  shall  be  employed  to  assess  and  to 
diagnose  the  damaged  structure  and  to  describe  the  extent  of 
damage  from  the  standpoint  of  energy  criterion.  In  addition,  the 
established  structural  power  can  also  be  a  material  character — 
istic  for  nonlinear  structures.  Such  characteristic  is  important 
for  identification  of  nonlinear  structure  as  well  as  damage. 

The  fundamental  purpose  of  structures  is  to  protect  or 
shield  the  occupants.  However,  while  structures  may  be  designed 
adequately  for  the  anticipated  excitation,  many  develop  serious 
structural  damage  due  to  unexpected  severe  loadings  or  deteriora¬ 
tion  from  exposure  to  elements.  As  a  consequence,  the  structure 
may  behave  nonl inear ly  during  excitation.  Furthermore,  for  most 
cases  the  nonlinear  behavior  occurs  even  for  undamaged  status. 

The  nonlinearity  must  be  characterized  in  order  to  predict  the  - 
structural  response. 

Currently,  the  available  techniques  and  theories  for  damage 
diagnosis  are  (i)  visual  observation  <ii)  volumetric,  such  as 
X-ray,  radiography  and  magnetic  field  methods,  (iii)  dynamic  me¬ 
thods,  such  as  the  acoustic  wave  and  modal  theories.  The  first 
two  methods,  visual  and  volumetric,  are  essentially  local,  thus 
not  applicable  to  hard-to-reach  locations  or  to  complex  struc¬ 
tures.  Hence,  they  cannot  be  adapted  to  establish  the  safety  and 
reliability  of  a  structure  for  future  excitation.  In  other  words, 
the  engineer  cannot  accurately  estimate  the  reliability  for  the 
damaged  structure  without  an  analytical  theory  to  support  a  relia 
bility  analysis.  The  dynamic  theories  are  based  upon  the  measure 
ment  of  structural  dynamics  characteristics  prior  to  and  after 
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the  dmge  [19],  or  upon  tht  recording  of  structural  response  dur¬ 
ing  a  known  strong  excitation  £13,16,17,213.  The  dynamic  theories 
are  known  to  have  uncertainties.  A  probabilistic  theory  is  being 
developed.  The  strong  excitation  diagnostic  theory  depends  on  an 
accurate  nonlinear  system  identification. 

The  material  nonlinearity  for  most  cases  is  characterized  by 
its  force  deformation  relationship.  For  example,  the  force-def or- 
mat ion  hysteretic  curves  for  elasto-plastic ,  work-hardening, 
strength  degradation,  or  softening  materials  display  different 
shapes  which  may  be  displacement  dependent.  If  the  force-deforma¬ 
tion  relationship  is  known  precisely,  it  is  possible  that  engin¬ 
eers  can  predict  the  structural  response  with  good  accuracy. 

The  problem  here  is  that  can  we  know  precisely  the  nonlinear 
force-def or mat i on  behavior?  For  static  case,  the  answer  probably 
is  positive,  since  the  static  assumption  excludes  the  frequency 
dependent  property.  The  material  then  may  behave  closely  as 
those  obtained  in  lab  tests.  However,  when  transient  analyses 
are  desired,  many  structures  may  or  may  not  follow  the  rules 
which  are  developed  in  accordance  with  the  lab  test.  For  such  a 
case,  the  difficulty  is  to  measure  the  material  restoring  force 
during  extreme  excitation  when  structure  behaves  nonlinearly. 

In  many  cases,  it  is  extremely  difficult  to  measure  the 
structural  hysteretic  restoring  force  due  to  the  environment 
restriction.  As  a  consequence,  the  predicted  structural  response, 
using  lab  developed  material  model,  doesn't  correlated  very  well 
with  the  measured  response,  such  as  acceleration.  The  present 
paper,  therefore,  aims  to  establish  an  alternate  approach  to 
characterize  the  nonlinear  structural  behavior.  The  established 
material  characteristic  can  be  employed  to  define  the 
nonlinearity  of  the  structure  during  a  dynamic  response,  to  model 
the  structural  system  behavior,  and  to  predict  the  damage  status 
of  the  structure. 
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#  .  a  rorwuiuon  ut  otruttur»i  rower  rmue;  i 


The  nonlinear  system  considered  herein  is  of  single -degree-of — 
■freedom  <SDF) ,  shown  in  Figure  9.1,  with  the  governing  differen¬ 
tial  equation  of  motion 


m  z  +  c  z  ♦  R(z>  *  fit). 


a 


Rcz)  ,  nl 


— yWvV*- — f 


Fig.  9.1  The  SDF  System 


(9.1> 


where  the  mass  <m)  and  the  damping  ratio  (c)  are  deterministic, 
fit)  is  the  random  excitation,  it,  z,  z>  are  respectively  the  dis¬ 
placement,  velocity  and  acceleration  responses  of  the  system,  Riz) 
is  the  displacement-related  nonlinear  system  hysteretic  force.  The 
hysteretic  restoring  force  is  characterized  by  the  materials, 
which  can  be  clasto-perf ect  plastic,  work-hardening,  or  stiffness/ 
strength  degradation  £13,17,213.  The  property  of  stiffness/ 
strength  degradation,  shown  in  Figure  9.2,  is  the  typical  behavior 
of  concrete  structures  subject  to  strong  excitation.  As  a  result, 
the  energy  dissipates  through  the  hysteretic  loop.  The  energy 
dissipation  (u)  is  thus  defined  ass 


Fig.  9.2  Typical  Behavior  of  System  Degradation 
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u  * 


I 


z  (t) 

o 


R<f)  df. 


<9. 2) 


which  represents  the  arcs  under  the  hysteretic  loop.  It  has  been 
deaonstrated  in  1133  that  the  amount  of  energy  dissipation  can  be 
used  as  an  indicator  for  the  damage  level,  which  is  associate  with 
the  diminishing  of  strength  and  stiffness  of  the  structure.  How¬ 
ever,  currently  the  amount  of  energy  dissipation  is  very  difficult 
to  measure,  especially  for  a  continuous  structure.  In  view  of 
this  problem,  we  introduce 


q(t)  *  du/at. 


(9.3) 


which  represents  the  energy  dissipation  rate.  Hence,  the  equation 
of  motion  <9.1)  becomes: 

m  z  *  c  z  +  q(t)/z  *  f <t).  (9.4) 

Since  q(t)  represents  the  rate  at  which  the  energy  is  dissipated 
through  the  hysteretic  loop,  it  can  be  viewed  as  a  structural 
power .  which  is  an  important  characteristic  of  a  structure.  If 
q(t)  can  be  expressed  in  terms  of  measured  responses,  damage  then 
can  be  identified  by  using  q(t)  as  a  damage  indicator.  Therefore, 
if  we  introduce  the  normalized  kinetic  energy  (T) ,  such  that 


T  *  z2/2, 

the  equation  of  motion  (9.4)  is  thus  formulated  as: 


(9.5) 


2c T  +  q (t )  *  f (t)  z  -  a  8T/8t.  (9. 6) 
The  formulation  (9.6)  can  be  applied  to  the  following  problems: 

(i)  For  an  undamaged  structure,  the  value  of  the  structural 
power  can  be  identified  by  computing  the  velocity  and  a  con 
trolled  excitation.  Such  process  constitutes  the  nonlinear 
system  identification. 
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<ii>  Fear  ■  dAMQtd  structure,  ths  value  of  the  structural  pow¬ 
er  will  dieinish  as  compared  to  that  of  the  undamaged  struc¬ 
ture.  In  order  to  assess  the  extent  o-f  damage,  the  value  of 
the  structural  power  of  a  damaged  structure  must  be  determin¬ 
ed.  It  can  be  done  by  using  the  process  as  discussed  in  (i). 

< i i i  >  I-f  the  velocity  response  is  monitored  during  field 
excitation,  such  as  earthquake,  the  value  of  the  structural 
power  can  be  estimated  during  the  course  of  the  excitation. 
With  the  estimated  value  of  q(t>,  the  damage  status  at  the 
end  of  excitation  can  be  readily  enumerated. 


9.2  Pi scussi on  of  Structural  Power t 

The  structural  power,  as  defined  in  Sec.  2.0,  is  frequency  de 
pendent.  To  demonstrate  the  characteristic  of  frequency  dependen¬ 
cy,  the  Q-hysteresis  model,  one  of  the  proposed  models  for  the  be¬ 
havior  of  stiff ness/strength  degradation,  is  adopted  here  to  simu¬ 
late  the  structural  response.  The  excitation  illustrated  first  is 
a  narrow  band  cosine  wave  random  excitation: 

f<t>  «  A  cos  <uft+4),  C9.7) 

where  A  is  the  amplitude,  the  phase  angle  I  is  a  random  variable 
with  uniform  distribution  between  0  and  2s.  When  the  system,  as 
in  Figure  1,  with  the  restoring  force  characteristic  shown  in  Fig¬ 
ure  9.3,  is  excited  by  the  loading  given  by  (9.7),  the  energy  dis¬ 
sipation  is  computed,  for  three  different  excitation  frequencies 
—  123,  61.5,  30.75  rad/sec.  Figure  9.4.  The  figure  shows  that 
the  structural  power,  the  slope  of  the  energy  curve,  is  constant 
for  a  specific  narrow  band  random  excitation. 


47 


Fig.  9.3  The  Q-hysteresis 


Fig. 9. 4  Energy  Dissipation  for 
Narrow— band  Excitation 


To  characterize  the  structural  power  of  a  systee  subjected  to 
Mide  band  random  excitation,  the  excitation  Mill  be  expressed  as: 

N 

f<t>  *  E  A.  cos(u  t+4. ),  (9.8) 

i«-N  1  11 

where  Ai t  i*-N,..,N,  are  constant,  are  the  frequencies  uniform¬ 
ly  distributed  between  the  frequency  spectrum  range  (-100,  100)Hz, 
and  the  phase  angles  are  random  variables  with  uniform  distribu¬ 
tion  between  (0,2x) .  With  the  same  structural  system,  the  normal¬ 
ized  energy  dissipation,  which  is  defined  as  the  energy  dissipa¬ 
tion  from  G-hysteresis  normalized  by  the  quarter  cycle  energy  dis¬ 
sipation  at  failure,  was  computed  for  50  different  excitations,  gi¬ 
ven  by  the  wide  band  excitation  (9.8).  The  results  are  shown  in 
Figure  9.5,  from  which,  the  structural  power  can  be  assumed  as 
constant  with  the  following  statistics: 

ECq(t)3  “  0.116  S/%  •  constant. 

Standard  deviation  of  q  *  0.011  j/s. 

Furthermore,  in  order  to  illustrate  the  profile  of  the  struc¬ 
tural  power  vs  the  frequency,  the  same  structural  system  is  excit¬ 
ed  by  the  excitation  as  given  in  (9.7),  with  frequency  ranged  from 
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0.  0.1  l.t  l.l  !.«  J.O 

TltC  (SCC) 

Fig. 5  Energy  Dissipation  Sequence  -for  wi de-band  Excitation 

0.05  to  50  Hz.  The  results,  shown  in  Figure  9.6,  illustrate  that 
the  structural  power  (q)  decreases  as  the  excitation  -frequency  in¬ 
creases.  The  aaxieua  value  of  q  occurs  when  the  excitation  fre 
quency  is  in  the  neighborhood  of  the  natural  frequency. 


Finally,  in  order  to  demonstrate  the  influence  of  the  ampli¬ 
tude  of  the  excitation  on  the  structural  power ,  the  same  structure 
system  is  excited  by  a  sequence  of  different  amplitude  loadings. 
Figure  9.7  plots  for  3  different  frequencies  <0.5,  1.0,  3.0  Hz) 
the  resulting  structural  power  vs  the  excitation  aaplitude,  normal¬ 
ized  by  the  yielding  strength.  It  readily  shows  in  Figure  9.7, 
as  excitation  amplitude  rises  the  value  of  the  structural 
power  also  increases. 


Forcing  Frequency  (r/c) 

Fig. 9. 6  Structural  Power  vs 
Excitation  frequency 


Excltotlon  Amplitude 

Fig. 9. 7  Structural  Power  vs 
Excitation  Amplitude 
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Th«  technique  of  using  structural  power  to  assess  the  nonlin¬ 
ear  structure  requires  the  measurement  of  the  structural  velocity 
response  which  can  be  obtained  by  integrating  the  acceleration  res¬ 
ponse.  Therefore,  it  can  be  applied  to  the  daaage  assesseent  con¬ 
veniently.  It  was  demonstrated  that,  generally,  the  structural  pow¬ 
er  possesses  a  constant  value  for  a  specific  excitation,  either 
narrow  or  wide  band.  The  characteristics  of  constant  value  can  be 
further  assumed  as  a  random  variable  if  the  excitation  is  random. 
Such  assumption  simplifies  the  random  analysis  for  structural  res¬ 
ponse.  The  application  of  using  structural  power  in  random  vibra¬ 
tion  analysis  is  discussed  in  £233. 
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